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Abstract: In order to deepen our understanding of the nature of the deconfinement 
phase transition for various gauge groups, we investigate SU(4) Yang-Mills theory in 2 + 1 
dimensions. We find that the transition is weakly first order. We perform extensive Monte 
Carlo simulations on lattices with temporal extent Nt = 3, 4, and 5, and spatial sizes up 
to N s = 20 Nf. We observe coexistence of confined and deconfined phases at the critical 
temperature, and finite-size scaling shows consistency with first order exponents. The 
continuum extrapolation of the latent heat yields L^/T^ = 0.188(17). 
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1. Motivation 

The confinement of quarks inside baryons and mesons is a feature of the strong interactions 
only at low temperatures. At sufficiently high temperatures hadrons melt and quarks and 
gluons form a plasma. Numerical simulations of QCD support the expectation that this 
process takes place in a smooth manner and the high- and low-temperature regimes are 
analytically connected through a crossover. As one increases or decreases the quark masses, 
the situation changes and, at some point, a finite temperature phase transition occurs. This 
phase transition is related to symmetries that are either badly broken or only approximate 
when the quarks have their physical masses. For massless quarks, chiral symmetry is 
exact: it is spontaneously broken at low temperatures and it gets restored in a chiral phase 
transition at finite temperature. In the opposite limit, when the quarks are heavy, they 
are only weakly coupled to the gluons. As the quark mass becomes much larger than the 
typical energy scale of the strong interactions, Aqcd ~ 250 MeV, the quarks decouple and 
gluons are the only relevant degrees of freedom. In the limit of infinitely heavy quarks, the 
global center symmetry of Yang-Mills theory is no longer explicitly broken by the quarks' 
triality and it becomes an exact symmetry of the theory. When the temperature is about 
300 MeV, the center symmetry breaks spontaneously at a phase transition [1-3]. Since 
quarks transform non-trivially under center transformations, the breaking of the center 
symmetry implies deconfinement. 

Yang-Mills theory provides a simplified framework in which the phenomenon of con- 
finement can be investigated without facing the more difficult numerical problems related 
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to dynamical fermions. Moreover, the pure gluon dynamics is quite rich and many inves- 
tigations have been performed to study various interesting features of Yang- Mills theory. 
For instance, string effects in the static quark potential have been observed numerically [4]. 
The study of topological objects and of their relevance for the mechanism of confinement 
is an active field of investigation [5]. A systematic study of SU(N) Yang-Mills theory for 
various N is a research topic that aims at understanding the way in which the large N limit 
is approached [6, 7]. This paper deals with another important characteristic of Yang-Mills 
theory, namely the order of the deconfinement phase transition. 

About 25 years ago, Svetitsky and Yaffe conjectured [8] that the critical behavior of 
a gauge theory at the deconfinement transition can be described by a scalar field theory 
with a symmetry corresponding to the center of the gauge group. In fact, if one integrates 
out the spatial components of the gluon field in {d + l)-dimensional Yang-Mills theory, 
one obtains an effective action for the scalar field represented by the Polyakov loop. The 
corresponding scalar field theory is defined in d dimensions and, in general, its action is very 
complicated. However, if the deconfinement phase transition happens to be second order, 
as one approaches the critical point, the correlation length diverges and universal critical 
behavior arises. Hence, the details of the complicated effective action become irrelevant: 
only the center symmetry and the dimensionality of space determine the universality class. 

Svetitsky and Yaffe's conjecture has been checked in many numerical simulations in 
Yang- Mills theory with various gauge groups, both in 2+1 and in 3+1 dimensions. In those 
cases in which the deconfinement phase transition is second order, the universality class 
has indeed turned out to be the one predicted by Svetitsky and Yaffe. In 3 + 1 dimensions, 
SU{2) Yang-Mills theory has a second order deconfinement phase transition [9-14] in the 
universality class of the 3-dimensional Ising model [15, 16]. However, in 3 + 1 dimensions 
no other pure gauge theory has been found to have a second order deconfinement phase 
transition [17-26]. In 2 + 1 dimensions, SU{2) Yang-Mills theory again has a second order 
deconfinement phase transition, now in the universality class of the 2-dimensional Ising 
model [27, 28]. Since in 2 + 1 dimensions fluctuations are stronger than in 3 + 1 dimensions, 
there are two more cases in which the deconfinement phase transition is second order. At 
its deconfinement phase transition, (2 + l)-dimensional SU(3) Yang-Mills theory shows 
the same critical behavior as the 2-dimensional 3-state Potts model [28-30]. The group 
Sp(2) has the same center Z(2) as SU(2), and the deconfinement phase transition of the 
corresponding (2 + l)-d Sp(2) Yang-Mills theory is again in the universality class of the 
2-dimensional Ising model [31]. For other gauge groups, the deconfinement phase transition 
of the corresponding pure gauge theory has turned out to be of first order [32, 33]. 

In this paper, we examine SU(4) Yang-Mills theory in 2 + 1 dimensions because it is 
one of the last remaining unsettled cases. The original study on coarse lattices indicated 
that the transition is second order [34]. The improved numerical results presented in [35] 
show that, on coarse lattices, the transition is weakly first order. For finer lattices, the 
deconfinement phase transition appeared to be second order, perhaps belonging to the 
universality class of the 2-dimensional 4-state Potts model. This is a particular case of 
the 2-d Z(4)-symmetric Ashkin- Teller model which has lines of critical points along which 
the universality class and the critical exponents change continuously. The authors of [35] 
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pointed out that it is difficult to obtain a definite answer unless one considers rather large 
volumes and they could not rule out a weak first order phase transition in the continuum 
limit. The need for large volumes and fine lattices was also emphasized in [33]. 

In this paper we present numerical evidence for a weak first order deconfinement phase 
transition in 51/(4) Yang- Mills theory in 2 + 1 dimensions. Interestingly, although there 
is an infinite set of different available universality classes, the system deconfines with non- 
universal behavior. The center symmetry does not play a role in determining the order 
of the deconfinement phase transition. Only when the transition is second order does 
the center symmetry determine the universality class. As we conjectured in [31], the 
order of the deconfinement phase transition is determined by the size of the group. In 
the low-temperature confined phase, the dynamics of Yang-Mills theory is governed by 
glueballs. The number of glueball states — i.e. the number of singlets in the tensor product 
decomposition of adjoint representations — is essentially independent of the gauge group. 
On the other hand, the dynamics of the high-temperature plasma phase is determined 
by deconfined gluons, whose number is given by the number of generators of the gauge 
group. If there is a large mismatch between the number of relevant degrees of freedom in 
the confined and the deconfined phases, the phase transition does not proceed smoothly 
as a second order transition. Instead an abrupt discontinuous first order transition takes 
place. This conjecture is supported by numerical simulations which show that the strength 
of the first order transition increases with the size of the gauge group. Further evidence 
was provided by studies of Yang-Mills theory with the exceptional gauge group G{2) [36]. 
The group G{2) is the smallest, simply connected group with a trivial center. Therefore, 
in G(2) Yang-Mills theory there is no symmetry argument that implies the presence of 
a finite temperature deconfinement phase transition. However, since G{2) (which has 14 
generators) has a rather large size, G(2) Yang-Mills theory has a first order deconfinement 
phase transition in 3 + 1 dimensions [37-39]. Various aspects concerning the problem of 
confinement in G{2) Yang- Mills theory have been investigated in [40,41]. 

The rest of the paper is organized as follows. In section 2, we describe the standard 
lattice formulation of Yang-Mills theory, the observables that we consider, and the finite- 
size scaling analysis used to determine the order of the phase transition. The numerical 
results are presented in section 3, followed by our conclusions. 

2. SU(A) Yang-Mills Theory on the Lattice 
2.1 The action and the observables 

We perform numerical simulations of SU(4) Yang-Mills theory on a periodic lattice in 2 + 1 
dimensions. We consider the standard Wilson plaquette action 

S M = -j E ReTr u n = -j E ReTr (U^U x+ ^ul + ^ul v ), (2.1) 

where the link parallel transporter matrices U x „ £ SU(4) are group elements in the fun- 
damental representation. All dimensionful quantities are expressed in units of the lattice 
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spacing. The bare dimensionful gauge coupling (3 is related to the usual gauge coupling g 
in the continuum by = 8/g 2 . The path integral measure and the partition function Z 
then take the form 

[VU = T[[ dU x ^ Z= fvUexp(-S[U]). (2.2) 

The Polyakov loop [1, 2] 

1 Nt 

<Z> S =-Tr(Pl[Ug M+1 ) (2.3) 
t=i 

is the trace of a path-ordered product of link variables along a loop wrapping around 
the periodic Euclidean time direction. Here N t = l/T is the extent of the lattice in 
the Euclidean time direction, which determines the temperature T in lattice units. The 
expectation value of the Polyakov loop is given by 

($) = \ / ^^E^^m^)' ( 2 - 4 ) 

X 

where N s is the extension of the lattice in the spatial directions. The Polyakov loop rep- 
resents a scalar field that transforms non-trivially under symmetry transformations in the 
center subgroup Z(4) of SU(4). Hence, a non- vanishing expectation value of the Polyakov 
loop indicates the spontaneous breakdown of the center symmetry and thus signals decon- 
finement. However, in a finite periodic volume spontaneous symmetry breaking — in the 
sense of a non-vanishing order parameter — cannot occur. Therefore, in the finite-size 
scaling analysis discussed below, we will consider the expectation value of the magnitude 
of the Polyakov loop ( | | ) . In a finite volume this quantity is always non- vanishing but it 
approaches zero when one takes the thermodynamic limit in the confined phase. Another 
quantity that is useful for distinguishing the confined from the deconfined phase is the 
probability distribution for the Polyakov loop, 

V^) = ^ J P^^^-^^^Tr^n^+O^exp^C/]). (2.5) 

In the confined phase £>(<£) has a single peak centered at $ = 0. In the deconfined phase, 
on the other hand, it has four degenerate maxima at $ = <l?o exp(i/c7r/2), where <3?o is a 
positive real number and k = 0, 1,2,3. When the deconfinement phase transition is first 
order, the confined and the deconfined phases coexist and can be distinguished by their 
different values of the Polyakov loop even at the phase transition. In that case, close to 
the phase transition one thus observes five maxima of the distribution p(<3?). The relative 
weight of the confined and deconfined peaks changes as one crosses the phase transition. On 
the other hand, when the deconfinement phase transition is second order, the high- and 
low-temperature phases become indistinguishable at criticality. The confined maximum 
becomes broader and broader as the critical temperature is approached from below and, 
at criticality, the width of the peak diverges. When the temperature is increased further, 
the four Z(4)-symmetric deconfined peaks emerge smoothly from the broad distribution of 
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the Polyakov loop. In the limit of very high temperatures, one can perform an analytic 
perturbative calculation of the effective potential for the Polyakov loop [42-46] . 

Other useful observables that characterize the deconfinement phase transition are the 
Polyakov loop susceptibility %, defined by 

x = N * (m 2 ) - mf) , (2.6) 

and the specific heat C given by 

For a first order transition it is interesting to also consider the latent heat L^. As we 
said above, in the case of a first order transition one can distinguish the confined from the 
deconfined phases even at the transition. One can then define the action densities s c and 
Sd for the confined and the deconfined phases, respectively. The latent heat is defined as 
the difference in the action density between the two phases at the critical temperature and 
it is given by 

L h = s d - s c . (2.8) 

The fluctuations in the action attain their maximum when the two phases have the same 
probability. It then follows that, in the thermodynamic limit, the maximum C max of the 
specific heat is given by 

^max = ^2^ L 2. (2.9) 

For a second order deconfinement phase transition the latent heat vanishes since the con- 
fined and deconfined phases become indistinguishable at the critical point. 

2.2 Finite-size scaling 

Away from a phase transition, the susceptibility of an extensive quantity scales with the 
volume. This scaling behavior changes as we approach a phase transition and the fluc- 
tuations become stronger. For a first order transition, the susceptibility of an extensive 
quantity increases with the square of the volume. This scaling behavior follows from the 
coexistence of the two phases. In fact, in general, an observable has different values in the 
two phases. Since the observable is extensive, its susceptibility scales with the square of 
the volume. In case of a second order phase transition, the susceptibility scales faster than 
the volume but — unlike for a first order transition — not as fast as the square of the 
volume. The exponent that characterizes the scaling behavior depends on the observable 
and on the universality class of the phase transition. 

The method we have used in the finite-size scaling analysis of our numerical data is 
the following. For a given lattice size N% x Nt, we perform a set of numerical simulations 
at various couplings f3 across the deconfinement phase transition. Using the Ferrenberg- 
Swendsen re- weighting technique [47,48] we determine a pseudo-critical coupling (3 c ,N s ,N t 
from the maximum of the Polyakov loop susceptibility. We then repeat this procedure for 
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various values of the spatial lattice size N s . For a first order phase transition the critical 
coupling depends on the spatial volume as 



N? 

f3c,N s ,N t = 0c,N t + a 0jp + • • • (2.10) 

where /? c ,iv t is the critical coupling in the limit of infinite spatial volume at fixed temporal 
extent Nt. 

Up to corrections to scaling, the data for the Polyakov loop susceptibility density 
x/Ng collected at different couplings (3 and for different lattice sizes N 2 x N t collapse onto 
a single universal curve once they are plotted as a function of the finite-size scaling variable 
x = (N s /N t ) 2 (f3/(3 Ct N s ,N t — !)• The corrections to scaling can then be easily measured at 
the maximum of the curve 

v max / v max \ at2 

1 ' ' +bo^ + --. (2.11) 



TV? V /oo X 2 

S \ S / oo s 

A similar formula holds for the maximum of the specific heat 



^max i i\t2 

1 Ll + c ^ + ... (2.12) 



3N*N t 4 h ' u iV| 



3. Discussion of the numerical results 



3.1 Simulation details 

We have performed simulations on lattices with Nt = 3, 4, and 5, and for spatial sizes 
N s as large as 20 N t . We have used a standard combination of heat-bath [49] and over- 
relaxation [50-53] algorithms to update the various SU(2) subgroups of SU(4) [54]. We 
have simulated with a ratio of over-relaxation to heat-bath updates of 4/1 and 1/1, and 
we find no significant difference between these two choices. For each set of /?, Nt, and N s , 
we have generated at least 10 5 configurations to be used for measurements. These runs 
are sufficiently long such that for the smaller physical volumes with e.g. N s /N t = 10, we 
see 0(50-100) tunneling events between the various bulk phases. For the larger volumes 
like N s /Nt = 20, this is reduced to O(10) tunnelings. We expect that this is a reasonable 
sampling of the various bulk phases. 

3.2 Monte Carlo histories and Polyakov loop distributions 

In Figure 1, we plot the Monte Carlo histories of the Polyakov loop and the plaquette 
expectation value, configuration by configuration, for a 60 2 x 3 lattice at = 20.40. The 
system spends a long time in a particular phase, characterized by the value of <E>, before 
it rapidly tunnels to a different phase, in which it again remains for a significant period 
of Monte Carlo time. In this particular run, we see the confined phase, in which <I> fluc- 
tuates around 0, and four deconfined phases, in which $ varies around the four values 
<3?o exp(i/c7r/2). The lower plot shows that the plaquette changes simultaneously with the 
Polyakov loop, between two similar but still distinguishable values. This suggests that the 
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system is close to the phase transition, and that the deconfinement transition is first order, 
with coexisting confined and deconfined phases. However, the small jump in the plaquette 
value indicates that the transition may well be rather weak. This simulation is quite typical 
in that we see about 10 tunneling events occur in this large volume. 

In Figure 2, we show the probability distributions of 
the complex- valued Polyakov loop for simulations on 36 2 x 3 
lattices at three different (3 values, ranging from low to 
high temperature. At low temperature, there is just the 
confined phase, while at high temperature there are four 
deconfined phases. At (3 = 20.26, the system is apparently 
quite close to the transition temperature and the five bulk 
phases coexist. Figure 3 shows the probability distributions 
of the plaquette value for the same simulations. We find 
only a single-peak distribution, which varies smoothly with 
j3. In the case of a normal-strength first order transition, 
close to the critical temperature one would expect to see 
two distinct peaks. In the present case, the discontinuity 
in the plaquette value is clearly visible in the Monte Carlo 
history — as shown in the bottom part of Figure 1 — but 
due to its small size, it does not stand out in the plot of 
the probability distribution. 

3.3 Polyakov loop susceptibility 

We use re-weighting of the various ensembles to determine 
the location (3 c ,N s ,N t of the peak of the Polyakov loop sus- 
ceptibility x- This is an accurate method and we show 
typical results in Figures 4 and 5. We use the bootstrap 
method to calculate both the error in % f° r each individual 
ensemble, and the error in f3 c ,N a ,N t extracted from the re- 
weighting of several ensembles. We list the finite-volume 
Table 1- The finite- volume cr hical couplings in Table 1 for the various simulations we 
critical couplings for various val- have performed. 

ues of N s and N t , with boot- If the deconfinement transition is of first order, the 

strap errors. infinite- volume critical coupling should be approached asymp- 

totically as (N t /N s ) 2 , as described in Equation (2.10). The 
data as shown in Figures 6, 7, and 8 display exactly this behavior. In determining the 
infinite- volume critical coupling for the various Nt values, there is some systematic error 
involved in the choice of the extrapolation range. One does not know how large N s has to 
be before linear behavior sets in. Let us first discuss the N t = 3 data. We start by fitting 
all of the data, then omit one at a time the data for the smaller values of N s . The results 
are listed in Table 2, with a statistical error included for each fit. For Nt = 3, we see that 
the fit does not improve as data is omitted, so the optimal fit includes all of them. For 
the final error estimate, we use the jackknife method applied to the optimal data set. We 
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Pc,N s ,N t 


3 


26 


20.251(8) 
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20.271(6) 
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20.288(7) 
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20.3230(25) 
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20.351(5) 
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25.966(14) 




36 


26.019(14) 
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31.64(4) 




40 


31.75(4) 




46 


31.72(4) 
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31.77(3) 
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31.91(4) 
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31.977(26) 




66 


32.012(24) 




100 


32.090(22) 
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quote a final result of (3 c ,N t = 20.414(5). The optimal fit and this extrapolated value are 
also shown in Figure 6. 
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32.22(8) 





Table 2: The infinite- volume critical couplings for 
the various values of N t and the quality of the linear 
extrapolations, where some of the smaller N s data are 
discarded. 



We apply exactly the same proce- 
dure to the Nt = 4 results. From Fig- 
ure 7, we see that a linear fit of all 
of the data gives (3 c ,N t = 26.251(13), 
the error being statistical only. How- 
ever, the quality of the fit is poor, with 
X 2 /d.o.f. = 17.4/6. Omitting the data 
for smaller N s , the quality of the fit 
does not improve, in fact it becomes 
worse, as listed in Table 2. The small 
N s data do not seem to be at fault, as 
significant fluctuations occur at larger 
N s . A quadratic fit gives (3 c ,N t = 26.26(3) 
and x 2 /d.o.f. = 17.2/5, so the data show 
almost no quadratic behavior and the 
extrapolated value agrees very well with 
the linear fit. Hence, the optimal choice 
seems to be to use all of the data. Per- 
forming a jackknife error analysis on 
this set, the final result we quote is (3 c ,N t = 
26.251(16). 

For both Nt = 3 and 4, a linear ex- 
trapolation appears valid for a fitting 
range (N t /N s ) 2 < 0.015. We expect 
that the same should hold for N t = 5. 
Looking at the data in Figure 8, two of 
the data lie outside this range. At suf- 
ficiently small N s , one does expect to 



see deviations from linearity. A linear fit of all of the data gives a very poor quality of fit. 
When the two smallest N s values are omitted, x 2 /d.o.f. improves significantly, from 26.1/6 
to 7.8/4, giving /3 c ,N t = 32.22(4), the error being statistical only. As listed in Table 2, 
the fit improves further if the next two small N s data are discarded. However, we believe 
that this is excessive and gives the impression of a more accurate extrapolation than is 
warranted. Using the linear regimes of Nt = 3 and 4 as a guide, we conclude that the 
optimal fit excludes only the two smallest N s data. For comparison, we also fit all of the 
data using a quadratic form, which gives P c ,N t = 32.24(6) and x 2 /d.o.f. = 14.9/5. The ex- 
trapolated value is completely consistent with that from a linear fit, although the quadratic 
fit is somewhat inferior. Using the jackknife method on the optimal set, we quote a final 
value for N t = 5 of (3 CtNt = 32.22(8). 
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We also investigate the finite-size scaling behavior of the Polyakov loop susceptibility. 
In Figures 9, 10, and 11, we plot the rescaled susceptibility x/^s as a function of the 
finite-size scaling variable, x = (N s /N t ) 2 (/3/(3 c ,N s ,N t — 1), for N t = 3,4, and 5. We are 
assuming here that the exponents are those of a first order transition. We see that, close 
to the critical temperature, the data do indeed collapse onto a single universal curve. This 
is further evidence of the first order nature of the deconfinement transition. 
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0.635(13) 
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1.091(5) 
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56 


0.739(3) 




60 


0.646(3) 
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0.5352(29) 




100 


0.2406(26) 



Table 3: The finite- volume maxima of the 
specific heat for various values of N s and 
N t , with bootstrap errors. 



3.4 Specific and latent heats 

Besides the Polyakov loop, we wish to use an- 
other thermodynamic quantity in order to deter- 
mine the order of the phase transition. We have 
attempted to measure the latent heat Lh by sim- 
ulating directly at the appropriate pseudo-critical 
coupling Pc,N a ,N t ■ Based on the value of the Polyakov 
loop, we have divided each ensemble into confined 
and deconfined configurations. We have then mea- 
sured the action in each phase and have deter- 
mined the discontinuity. Unfortunately, this method 
has some difficulties due to the somewhat arbi- 
trary cut in the value of the Polyakov loop used 
to distinguish confined from deconfined configu- 
rations. We have found that it is more accurate 
to measure (7 max ; the peak in the specific heat. 

In Table 3, we list the finite-volume maxima 
of the specific heat for N t = 3,4, and 5. The 
quoted errors for each ensemble are calculated 
using the bootstrap method. In Figures 12, 13, 
and 14 we plot the data for C max /(3iV s 2 iV t ) for 
N t = 3, 4, and 5 respectively, and for various val- 
ues of N s . It is no surprise that the data do in- 
deed extrapolate accurately in (Nt/N s ) , first 
order transition dictates. The infinite- volume val- 
ues C max / (SNg Nt) (oo) are quite small, and accu- 
rate data are needed in order to reach a reliable 
conclusion. Like the finite-volume critical cou- 
plings, one does not know a priori how large N s 
has to be before the linear regime in (Nt/N s ) 2 is 



reached. We find that omitting the smaller N s 
data does not improve the already very good linear fits, so all of the data are used in 
the final analysis. The results of the fits are presented in Table 4, where the errors are 
calculated by the jackknife method. 
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Using Equation (2.12), we convert the infinite- volume extrapolations into the latent 
heat Lh. We extrapolate the dimensionless quantity NfL^ = L^/T^ to the continuum 
linearly in 1/-A 2 , which describes the data very well. Our continuum determination is 
Lh/Tj? = 0.188(17), where the error is statistical, and the quality of the fit is x 2 /d.o.f. = 
0.90/1. With only three data points, it is not possible to estimate a systematic error by 
varying the choice of fitting range. The transition clearly becomes weaker on finer lattices, 
but its first order nature persists in the continuum limit. 



4. Conclusions 



Our results show that (2 + l)-d SU(4) Yang-Mills theory has a first order deconfinement 
phase transition. Large and fine lattices were important for reaching this result. In the 
various extrapolations and universal curves, our assumption of using first order exponents 
is well confirmed by the data. The first order nature of the transition is further supported 
by the observed coexistence of the confined and deconfined phases for all three temporal 
extensions N t = 3, 4, and 5 that we have considered. It is much harder to independently 
determine the critical exponents than to show consistency with an expected set. Since the 
2-d Z(4)-symmetric Ashkin- Teller model has continuously varying critical exponents, the 
challenge is particularly large. In [35], by including logarithmic corrections to scaling, the 
data suggested a second order deconfinement transition belonging to the universality class 
of the 2-d 4-state Potts model. The numerical data presented here suggest that this is 
not the case. The determination of a non-zero latent heat in the limit of vanishing lattice 
spacing shows that the deconfinement phase transition does not weaken to second order 
but stays first order in the continuum limit. 

It is surprising that so few gauge theories 
realize the Svetitsky-Yaffe scenario which only 
applies when the deconfinement phase transi- 
tion is second order. This may be particularly 
surprising in the present case, in which the in- 
Table 4: The extrapolated infinite-volume finite set of different universality classes of the 
maxima of the specific heat, and the quality 2-d Ashkin-Teller model would be available. In 
of the linear fits. Errors are calculated using 3+1 dimensions, SU{2) = Sp(l) Yang-Mills 
the jackknifc method. theory is the only pure gauge theory with a 

second order deconfinement phase transition. 
In 2 + 1 dimensions, the transition is second order only for SU{2), SU(3), and Sp{2) Yang- 
Mills theory. Even though all symplectic groups Sp{N) have the same center Z(2), the 
transition becomes first order as N increases, in both 3 + 1 and 2 + 1 dimensions. Inter- 
estingly, the 3-d Z(A)-symmetric spin model belongs to the universality class of the U (1)- 
symmetric XY model for N > 5, i.e. the symmetry is enhanced at the critical point [55]. 
However, the corresponding (3 + l)-d SU(N) gauge theory is unaffected by this peculiar 
critical behavior because its deconfinement transition is first order of a strength increasing 
with N. Indeed, as we conjectured in [31], not the center but the size of the gauge group 
determines the order of the deconfinement transition. In 3 + 1 dimensions all Yang-Mills 



N t 


C max /(3A s 2 A t )M 


rVd.o.f. 


3 


5.66(22) x 10~ 5 


5.8/8 


4 


5.5(5) x 10" 6 


0.65/6 


5 


1.19(11) x 10~ 6 


0.42/6 
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theories whose gauge group has more than three generators have first order transitions. 
As we now know, in 2 + 1 dimensions only the Yang-Mills theories whose gauge group 
has at most ten generators (namely SU(2), SU(3), and Sp(2), which has ten generators), 
have a second order deconfinement phase transition. An interesting case for future study 
is (2 + l)-d G{2) Yang-Mills theory. Since the exceptional group G(2) has a trivial center, 
there is no symmetry reason for a deconfinement phase transition and there may hence 
just be a crossover. In the absence of a non-trivial center, a second order phase transition 
can be ruled out on theoretical grounds, because it would require unnatural fine-tuning of 
some parameter. If the 14 deconfined G{2) gluons at high temperature cannot smoothly 
crossover to the low-temperature regime governed by a small number of glueball states, 
(2 + l)-d G(2) Yang-Mills should have a first order deconfinement phase transition. Since 
the 15 SU(4) gluons behave in this way, based on the size of the gauge group, one may 
expect the same for G{2). 
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Figure 1: The Monte Carlo histories of the Polyakov loop and the plaquette on a 60 2 x 3 lattice 
at j3 = 20.40, close to the deconfinement transition. The system tunnels between the confined and 
four deconfincd phases, with the plaquette value tracking the change in the Polyakov loop. 
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Figure 2: The probability distributions of the Polyakov loop on 36 2 x 3 lattices at j3 = 20.0 (top), 
20.26 (middle) and 20.5 (bottom). At low temperature, there is a single confined phase. Close to 
the critical temperature, we observe coexistence of the confined with the four deconfined phases. 
At high temperature, there are only the four deconfined phases. Because of the finite length of the 
simulation, the four deconfined phases are not equally sampled. 
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Figure 3: The plaquette value distributions for the 36 2 x 3 simulations shown in Figure 2. Close 
to the critical temperature, at = 20.26, we do not find two well-separated peaks, as one would 
expect for a normal-strength first order phase transition. 
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Figure 4: The Polyakov loop susceptibility x as a function of j3 for 42 2 x 3 lattices, using re- 
weighting of the combined ensembles. The critical coupling, where \ has a maximum, is determined 
to be (3 C} N s ,N t = 20.351(5). 
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Figure 6: The critical couplings f3 c ,N s ,N t for N t = 3 and a range of N s , extracted from the peak 
of the Polyakov loop susceptibility \- A linear extrapolation in (N t /N s ) 2 describes the data well, 
as expected for a first order transition. The infinite-volume extrapolated value is (3 c .N t — 20.414(5) 
(jackknife error) and the quality of the fit is x 2 /d.o.f. = 6.9/6. 
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Figure 7: The same as in Figure 6, this time for Nt = 4. The extrapolated value is {3 c ,N t — 
26.251(16) and the quality of the fit is x 2 /d.o.f. = 17.4/6. 
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Figure 8: The same as in Figures 6 and 7, this time for N t = 5. For N t — 3 and 4, the data 
show linear behavior for (N t /N s ) 2 < 0.015. Hence, the two smallest N s values are excluded from 
the optimal fit. The extrapolated value is (3 c ,N t — 32.22(8) and the quality of the fit is x 2 /d.o.f. = 
7.8/4. 
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Figure 9: The rescaled Polyakov loop susceptibility x/^s versus the finite-size scaling variable 
L 2 (P I 'f3 c ,N B ,N t ~ 1) f° r Nt — 3, where L — N s /N t . The data fall onto a universal curve, consistent 
with first order exponents. 
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Figure 10: The same as in Figure 9, this time for N t = 4. Again, the data are consistent with 
first order exponents. 
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Figure 11: The same as in Figures 9 and 10, this time for N t = 5. Again, the data are consistent 
with first order exponents. 
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(N t /N s ) 2 

Figure 12: The peak of the rescaled specific heat capacity C max /(3N^Nt) for Nt = 3 and 
a range of N s . A linear fit describes the data very well and gives the extrapolated value 
C ma 7(37V s 2 iV t )(oo) = 5.66(22) x 1CT 5 . The quality of the fit is x 2 / d -°- f - = 5.8/8. 
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Figure 13: The same as in Figure 12, this time for N t — 4. Again, we linearly extrapolate to 
obtain C max /(3N*N t )(oo) = 5.5(5) x 1(T 6 and the quality of the fit is x 2 /d.o.f. = 0.65/6. 
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Figure 14: The same as in Figures 12 and 13, this time for Nt = 5. A linear fit of all the data 
gives C max /(3iV s 2 iVt)(oo) = 1.19(11) x lfr 6 and the quality of the fit is % 2 /d.o.f. = 0.42/6. 
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Figure 15: The latent heat L h , taken from the N s -> oo extrapolation of C max / (3N^N t ). The 
data are extrapolated linearly in 1/Nf. The continuum value is NfLh = 0.188(17) and the quality 
of the fit is xVd.o.f. = 0.90/1. 
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